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I Abstract 

J^ . The propagation of Dyakonov surface waves (DSWs) at the planar interface 

f~^ I between an isotropic material and a linear electro-optic birefringent material 

OO ' can be dynamically controlled using the Pockels effect. The range of directions 

. ■ for DSW propagation has been previously found to be rather narrow. By careful 

^T ! choice of various parameters, this range of directions can be increased by more 

QQ ' than an order of magnitude. 
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1 Introduction 

Electromagnetic surface waves at the planar interface of two different materials travel 
along a direction parallel to the interfacial plane with decaying amplitudes perpen- 
dicular to the direction of propagation in both materials. Research on surface-wave 
propagation can be traced back to about a century ago to Zenneck pLj and continues 
to be of technoscientific importance [2] . 

In recent years, Dyakonov surface waves (DSW) have attracted much attention [3J. 
These waves can exist at the planar interface of an isotropic material and a uniaxial 
material, both dielectric, as was shown first by D'yakonov [3] and then enlarged upon 
by Averkiev and Dyakonov [5J. The isotropic material is not metallic. Later, several 
researchers considered DSW propagation on increasingly complex bimaterial interfaces 
such as biaxial-isotropic, uniaxial- uniaxial, and biaxial-biaxial materials — see, e.g., 
Refs. 6-10. 

DSW propagation is assured when a dispersion equation is satisfied. This equation 
arises from having to satisfy the frequency- domain Maxwell postulates in both mate- 
rials as well as the boundary conditions at the interface. Therefore, DSW propagation 
depends highly on the crystallographic symmetries of the two materials, at least one 
of which must be birefringent. Calculations for many pairs of materials have shown 
that conditions for DSW propagation are favored only for a narrow and restricted 
angular range of propagation directions in the interface [51 [3, dUl [HI [121 [IS]- The 
width of this angular existence domain (AED) is denoted here by Aij). With the ex- 
ception of some interfaces involving hypothetical metamaterials [12] , Aip for interfaces 
of non-magnetic materials are small — in most cases, only a fraction of a degree[3]. 
This makes experimental observation of DSW propagation very difficult. One of the 
largest values of Aip has been calculated to be ~ 4° [3j ; however, in that case one of 
the materials is calomel, a chalky and dusty material that is impractical to use. 

The potential real-life applications envisaged for the use of DSW are numerous [3], 
including devices for integrated optics, chemical and biological surface sensing, etc. 
To be able to experimentally launch a DSW, to experimentally detect a DSW, and to 
possibly control the DSW characteristics, a sufficiently wide AED appears desirable. 
That desire motivated this communication. 



In a predecessor paper [M], we looked at the control of DSW propagation at the 
interface of an electro-optic material, potassium niobate, and an isotropic material 
through application of an external dc electric field. Now, we examine optimizing the 
width Alp of the AED by an exploration of the effects of: the crystal orientation, 
the direction of applied dc electric field, and the strength of the applied dc electric 



field. Considering the planar interface of an isotropic and a birefringent material, 
both dielectric, we show in Sec. [2] that Aip can be increased by a large factor — if 
the birefringent material is susceptible to the Pockels effect — by suitable choices of 
parameters such as the magnitude and orientation of an applied dc electric field, the 
crystallographic orientation of the birefringent material in relation to the interface, 
and the refractive index of the isotropic material. 



2 Preliminaries and Numerical Results 

We considered the following problem p3]. The half-space z < is filled with a 
homogeneous, isotropic, dielectric material with an optical refractive index denoted 
by Ug. The half-space z > is filled with a homogeneous, linear, electro-optic material, 
whose optical relative permittivity matrix is stated as 

Cre/ = •S'^ (tp) ■ Ryix) ' ^PE " Ry{x) ' Sz ("V') • (1) 

Incorporating the Pockels effect due to an arbitrarily oriented but uniform dc electric 
field E_ '^, the matrix epE is given by 
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correct to the first order in \E '^l, where 
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are the principal relative permittivity scalars in the optical regime, whereas rjx 
(with 1 < J < 6 and 1 < -ft' < 3) are the electro-optic coefficients. The electro- 
optic material can be isotropic, uniaxial, or biaxial, depending on the relative values 
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of £-[ , €2 ; and 63 . Furthermore, this material may belong to one of 20 crystallo- 
graphic classes of point group symmetry, in accordance with the relative values of the 
coefficients tjk- 



The rotation matrix 
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'S'z(^) = sin Ip cos Ip (5) 



in Eq. ([T]) denotes a rotation about the z axis by an angle ip G [0, 27r). The matrix 

— sin X cos X 
Ryix) = I -1 I (6) 

cos X sin X 

involves the angle x ^ [0, 7r/2] with respect to the x axis in the xz plane, and combines 
a rotation as well as an inversion. The angles ip and x delineate the orientation of 
the electro-optic material in the laboratory coordinate system, the full transformation 
from laboratory coordinates {x, y, z) to those used conventionally for electro-optic 
materials (1,2,3) being illustrated by Nelatury et al. [1] 



The wave vectors on both sides of the interface were taken to lie wholly in the xz 
plane. The field representations based on the Maxwell curl postulates in the two half- 
spaces, the boundary conditions at the interface ^ = 0, and the dispersion equation 
for DSW propagation are available elsewhere [H]. For a simple case like that of the 
isotropic-uniaxial interface, obtaining a closed-form expression for Aip is relatively 
easy. However, for complex cases such as the one we are treating here, that goal is 
perhaps impossible. Recourse, then, has to be made to numerical methods in order 
to mount an extensive search for the proper choice of parameters that allow DSW 
propagation. 

A search was carried to maximize ^ip in relation to the magnitude and the direction 
of the dc electric field, the refractive index of the substrate, and the tilt angle x- 
Potassium niobate was chosen as the electro-optic material, because the magnitudes 
of its electro-optic coefficients are very large. 

Table 1 shows x ^-nd n^ values in a neighborhood wherein Aip is a maximum. 
Figure [1] shows the variation of A^ with n^ for x = 58.06°, 58.08° and 58.099°. The 
dc electric field is E^'^ = z 1.2 x 10^ V m~^. We found that, around x = 58°, even a 
small change in x produces a large change in A^. The last entry in Table 1 shows 
Alp ^ 1.4° which is 20 times larger than the result previously reported by us 



The search for a wider AED was also attempted choosing Ep'''^ parallel to —z, 
X cos Ip + z sin Ip and —x sin ip + z cos ip for the same choices of n^ and x values. But 
we could get only Aip of order 0.02°. However, we wish to continue the search hoping 
for wider AED. 



X 


Us 


At/- 


58.05° 


2.2805 


0.7312° 


58.06° 


2.2804 


0.7805° 


58.07° 


2.2804 


0.8179° 


58.08° 


2.2803 


0.9761° 


58.09° 


2.2802 


1.0302° 


58.099° 


2.2801 


1.4° 



Table 1: Maximum ^ip for some choice of parameters. 



3 Concluding Remarks 



In summary, Dyakonov surface waves propagate along interfaces involving birefringent 
materials over a narrow angular range A-?/; of the orientation angle ip. Having a larger 
A%Ij is helpful for easier generation and detection of DSWs. Several attempts have been 
made to widen this angular existence domain. By considering specific directions of an 
applied dc electric field and specific values of x, we have shown that exploitation of the 
linear electro-optic (Pockels) effect has the potential to widen the AED. The search 
for higher values of A-?/; for an arbitrarily oriented dc electric field is still open. Also, 
one might look into ranges of x values for specific values of ip and various magnitudes 
and directions of the applied dc electric field. 
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Figure 1: Variation of A-?/; versus n^ for different x values 
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